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We search for signals of Weak Annihilation in inclusive semileptonic D decays. We consider both 
the widths and the lepton energy moments, which are quite sensitive probes. Our analysis of Cleo 
data shows no clear evidence of Weak Annihilation, and allows us to put bounds on their relevance 
in charmless B semileptonic decays. 



I. INTRODUCTION 

The value of \Vub\ preferred by current global analyses of CKM data is about 15% smaller than the one extracted 
from inclusive charmless semileptonic B decays PJ. Though not very significant, the discrepancy has prompted a 
reexamination of the sources of theoretical uncertainty in the inclusive determination ,2^ . Weak Annihilation (WA) 
contributions are generally considered an important source of uncertainty in the Operator Product Expansion (OPE) 
that describes the inclusive B decays [3], and affect especially the high and lepton endpoint analyses. They appear 
in the OPE as I/ti^ corrections involving the matrix elements of dimension-6 four-quark operators, and affect both 
the total B — )• X^iv decay rate and the charged lepton energy spectrum [H |3] . 

Early estimates of the relevant matrix elements were derived in the framework of QCD sum rules They were 
also computed on the lattice in the static heavy quark limit [7] and with propagating heavy quarks [8]. However, 
it was soon realized [HI HO] that the light flavour component of the WA operators can differ from the light valence 
quark of the _B-meson leading to the so-called non-valence WA contributions which are extremely difficult to study 
non-perturbatively. In addition, contrary to valence WA contributions, the latter cannot be constrained by comparing 
charged and neutral B meson semileptonic decay rates and a common prejudice that the former should dominate may 
be unfounded. 

It was also noted in ref . pUHT^ that the WA matrix elements that enter B — ^ Xu^v decay can be constrained via the 
semileptonic decays of D and mesons, using heavy quark symmetry. Several authors have attempted to extract 
information on WA contributions from the measured total semileptonic rates of D^'^ and Ds, most recently in [13H15) . 
For instance, one may attribute the observed differences in and semileptonic widths [TB] 

T{D+ Xe+v)/T{D° Xe+iy) = 0.985(28) , T{Df Xe+iy)/T{D" Xe+v) = 0.828(57) (1) 

to the valence spectator quark WA contributions in Dg decays, since they are Cabibbo suppressed in the D'^ case and 
completely absent in D'^ decays. However, additional contributions to this difference arise from SU(i) breaking in the 
matrix elements of all dimension 5 and 6 operators, that contribute significantly to the total rates [TT]. Moreover, the 
slow convergence of the l/rric expansion, the generally large perturbative corrections to the semileptonic charm width, 
and the strong dependence on the charm mass may obscure the determination of the non- valence WA component 
from the widths ^. 

Our strategy in this paper is to consider also the moments of the lepton energy spectra in the OPE and compare 
them with recent experimental results from Cleo }16| . Not only are the moments free of the strong dependence on 
the charm quark mass and its associated uncertainty, but their perturbative and non-perturbative corrections tend to 
cancel as well. Moreover, since WA is expected to dominate the spectrum endpoint, leptonic moments might be more 
sensitive to WA contributions than the total rate. A study of moments in D semileptonic decays is also instrumental 
for a critical reassessment of the OPE in charm decays many years after jl7j , in view of the recent experimental results 
and of the successful application to B semileptonic decays. 

In the next Section, we present the experimental results on the leptonic spectra and compute from them the first 
leptonic moments. The OPE calculation of the moments is presented in Section HI, while we derive and discuss our 
results in Section IV. We conclude with a brief summary. 
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Figure 1: Normalized charm semileptonic decay spectra [16], extrapolated to pe = and boosted to the reference frame of the 
decaying particle. 



II. EXPERIMENTAL DATA 



Recently the Cleo Collaboration has measured the electron spectra of inclusive semileptonic and also Ds 

decays with a lower cut on the electron momentum in the lab frame of Pe > 0.2 GeV [TQ. They extract the total 
decay rates by extrapolating the spectra over the remaining phase-space, using a theoretically modeled sum over 
known exclusive modes. Unfortunately, they do not provide higher leptonic energy moments. 

Cleo give the electron energy spectra in the laboratory frame. Since the decaying D^'^ mesons are pair-produced 
with center-of-mass energy Ecm — 3.774 GeV, they are boosted in this frame with f3 — 0.14, 0.15, respectively. In the 
case of Dg mesons, the situation is more complicated since roughly half of them are produced in the primary vertex 
at Ecm = 4.170 GeV associated with D* mesons {(3 = 0.2). The other half then conies from radiative and hadronic 
D* decays. 

We extract the moments without cuts by first extrapolating the measured spectra towards Pe = 0. We use a 
procedure similar to the one employed by Cleo themselves, but instead of using a sum over exclusive modes with 
model form factors, we use the theoretical form of the spectrum in the « region, where the OPE is expected 
to provide a satisfactory description. In particular, for — > the spectrum must die at least like Pg. Thus we fit 
the first four measured bins to dT/dx ~ ax^{l + bx){l — x), with x — 2E^/mc. In the case of Z?s, due to the larger 
statistical uncertainties, we set 6 = and fit only to the first two measured bins. In the fits we take the reported 
systematic errors to be fully correlated between different bins. Using the fitted formulas we have computed the rates 
and obtained the following branching fractions 

B{D^ ^ Xev) = 0.064(1) , B{D+ Xev) ^ 0.161(3) , B{Ds Xeu) = 0.065(4) , (2) 

fully compatible with Cleo reported results. We then compute in the lab frame the first two leptonic energy moments 
normalized to the total rate, again assuming the systematic errors to be fully correlated among different bins and 
obtain 

(£;e)f, - 0.465(3)GeV , {El)f°, = 0.248(2)GeV2 , (3) 

(£;e)^.fc = 0.459(l)GeV , {El)f^, - 0.242(l)GeV2 , (4) 

{E,)f^, = 0.466(12)GeV , {El)f:, = 0.254(13)GeV2 . (5) 

The dependence of these values upon our extrapolation ansatz is expected to be smaller than the corresponding effect in 
the total rate, since the contribution of the pe ^ region to the n-th (unnormalized) moment is (2pe/m£))" suppressed 
compared to its contribution to the total rate. Thanks to the correlations between numerators and denominators, the 
total errors in ^ are significantly smaller than in ([2|. 

In order to compare these values to theoretical predictions in the D meson frame one still needs to take into account 
the boost factors. The energy of the final state electron in the lab frame is given by E'^ = ^E^(\ — f3cos9), where Eg, 
is the electron energy in the D meson frame and is the angle between the momentum of the electron and the one 
of the D meson in the lab frame. We thus obtain that {E'^) = {E^)i^^ = 7 {E^) with 7 = 1.009, 1.012 for Cleo's 
mesons, respectively. Similarly, for the 49% of I?s's which come from the primary vertex, 7 ~ 1.02. The remaining 
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51% receive an additional boost in an arbitrary direction since they originate from D* decays D* — )• Dg^ (94%) 
or £)* — DgiT (6%). In view of the precision in ([s]), the effect of this additional boost is always negligible. The 
directional averaging for the second electron energy moment {El) results in = (7^(1 + 1^'^ 1^)) {E'^) which again 

can be readily computed for all cases. The extrapolated rest-frame spectra are displayed in Fig. [l] Our results for 
the moments in the D mesons rest frame are 

{Ei)elp = 0.459(3)GeV , {E])f°^ = 0.240(2)GeV^ , ((S, - {Ee)f)f°j, - 0.029(2)GeV2 , 

= 0.455(l)GeV , {E])f^l = 0.236(l)GeV2 , ((S, - {Ee)^)^ = 0.029(l)GeV2 , (6) 



/ exp V / ' \ t / exp 

I exp 

^2 



= 0.456(ll)GcV , {E])f' = 0.239(12)GeV2 , ((S, - {E.yf)^:, = 0.031(12)GcV^ 



Here we have also listed the variances {<Te) of the spectra. Within the stated uncertainties there is no sign of a 
difference between the moments of Dg and D"'*. This is at odds with what one would naively expect from eq.([l]). 
However, as mentioned in the Introduction, the concurring contribution of SU(3) violation in the matrix elements of 
higher dimensional operators might provide a partial explanation. 

One can estimate the order of magnitude of SU{3) breaking corrections by comparing the hyperfine splitting of 
D^'^ and Dg mesons, A^^ = 3{m'jj, — m|,^)/4 (which is related to the OPE parameter ^q): 

A^-^+ = 0.409(l)GeV^ , A% = 0.413(l)GeV2 , = 0.440(2)GeV^ . (7) 

We see that even isospin violation between 1)+^° mesons is manifest, at the expected 1% level. On the other hand 
SU{3) violation is sizable and of the order 10%. Lattice QCD studies of fo and /d, il8] suggest that even 20% 
violation can be expected in certain quantities. Since the 1/to^ and 1/m^ contributions to the total semileptonic rate 
can be as large as 50% of the leading order estimate (depending on the charm quark scheme and scale), and since 
one can expect SU{3) violation of similar size in the matrix elements of the relevant power suppressed operators, an 
0(10%) SU{3) breaking in the widths is not unlikely. In order to reliably extract possible WA contributions from the 
measured total semileptonic rates of D and Dg mesons [TB], one would therefore need to estimate the size of SU{3) 
violation in all these matrix elements. In the case of normalized moments, on the other hand, some of the leading 
power corrections cancel out and one might be more directly sensitive to WA contributions. 

III. LEPTONIC MOMENTS IN THE OPE 

The perturbative corrections to moments of some kinematic distributions in semileptonic b decays are now known 
through 0(al) [5^. In the case of semileptonic charm decays, however, only the 0{as) ^25^ and 0{f3Qal) corrections are 
readily available using |21j . For what concerns the power corrections, explicit expressions for the leptonic spectrum 
at 0(l/m^) can be found in the last paper of [S], while we have computed the 0(l/m^) contributions from the 
form factors given in [22]. We neglect 0{l/m^) corrections [23] and 0{as/m1) corrections, for which only the n"^ 
contribution to the rate is known. 

We calculate the semileptonic rate and the total moments from the lepton energy spectra in terms of a; = 2iJe/TOc, 
where E^ is the electron energy in the D meson inertial frame. For simplicity, we only consider total leptonic moments, 
without lower cuts on the lepton energy, and neglect the lepton mass. We define 



Jo 



dx 1927r3 



TT TT'^ mi mi 



TO^ m 

,2 ^3 



(8) 



where r = m'^/m'^, ag = agirric), and /z^ g, j^g are the D meson matrix elements of the dimension 5 and 6 local 

operators appearing in the OPE. B^}^ is the WA contribution to the n-moment: as we will explain in a moment, it is 
related to the matrix elements of four-quark operators. In addition, there are also Cabibbo suppressed contributions, 
which can be included in the analysis by using the above formula in the limit r — > and by replacing Vcg with Vcd and 
similarly -Bw2'* with B^}'"^. They contribute to the total rate at the level of 5%, but their effect is highly suppressed 
in the normalized moments, with the possible exception of WA contributions. 

The tree-level spectrum in charm semileptonic decays is softer than in the analogous b decays and peaks at a; w 0.6, 
a feature qualitatively evident already in Fig. 1. The lowest order expressions for the rate and the first two moments 
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read 

/^°^ = l-8r + 8r^ -r* - 12r'^ log r , 

/^i) = [2r^ - iSr"* + eOr^ - 20r^ - 60r^ log r - 30r + 3]/5 , 

^ [_^6 ^ _ 3Q^4 ^ gQ^3 _ 35^2 _ gQ^2 log^ _ 24r + 2]/5 . 

The associated 0{as) corrections are known in closed form only for the total rate. We employ accurate numerical 
interpolations valid in the range 0.04 < y/r < 0.2 in the pole mass scheme 

/f ^ (r) = 2.86^^ - 3.84r log r - 2.47 , 
/}i)(r) = 2.11^^- 2.43rlogr- 1.56, 
/f ) (r) = 1.66V?- - l-71r log r - 1.09 . 

Similarly, while there is an analytic result for 72''^ (r) as an expansion in r |20j , we employ only the BLM approximation 
[21j in the form of numerical interpolations in the pole mass scheme: 

f^°\r) = /3o[8.16Vf - 1.21rlogr- 3.38] , 
/f^(r) = /3o[6.15Vf-0.35rlogr-2.24], 
(r) = ^o[4.99Vf + 0.1 Ir log r- 1.64], 

where /3o is the QCD beta function /?□ = 11 — 2n//3 and n/ = 3 is the number of light flavors. The above expressions 
can be translated to any other scheme (and associated IR scale) using the condition that the decay rate be scheme 
(and scale) independent at each perturbative order together with the known scale dependence of the heavy quark 
mass and power suppressed non-perturbative parameters. For the leading power corrections we get 

fi'\r) = (r)/2 , /(i)(r) = , /(^)(r) = ^/f (r) , 

= M - 2(1 - .)^ /«(.) = lf^'\r) - '(1 - r)^ fi!\r) = A/f (,) - ^(1 - r)^ 

/IS W = -f'J\r) , /W(r) = ^(1 - r)^ /f )(r) = ^/(^^(r) + ^(1 - 0^ 

fS\r) - ? + 0(r) + 81og ^ , /^(r) = ^ + 0(r) + 81og ^ , 
/l3''(r) = f + 0(r)+ 8 log 

where /iwA is the MS renormalization scale associated to the mixing of the Darwin and WA operators, an 0{a^) effect 

[261127]. The MwA dependence cancels against the implicit scale dependence of _B4rA '''(MwA)- A change of /iwA therefore 
shifts part of the Darwin operator contribution into the WA contribution, with important consequences for the error 
analysis. We find that for values of /iwA below IGeV the size of the coefficient in the width is comparable to that 
of other power corrections. 

WA contributions can be identified with the matrix elements of the dimension-6 four quark operators 

0?' = Q7m(1- 75)997^(1 -75)Q and Of'^Q{l~j,)qq{l^j,)Q, 
where Q is the heavy quark. The matrix elements that enter the charm meson decay rates are defined as 

Bl. = 2^ {D\ 01 - 0\ \D) = \niu!l{Bl^^ - J , (9) 

where ^ parameterize the deviation from the factorization approximation: Bwa vanishes in the limit of factorization. 
As it was recognized long ago [4j, WA is localized at the endpoint of the lepton energy spectrum, and can be 

approximately expressed by a delta function at the partonic endpoint. In this case, we would expect B^l'"^ — Swa, 
for all n, up to small 0{r) effects. In fact, gluon bremsstrahlung and hadronization effects are expected to smear 
the WA contribution over a region of electron energy around the partonic threshold (to^ — m1)/2mc corresponding 
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Figure 2; Different distributions of WA as a function of j: = 2Ee/mc, for rric = 1.4GeV. The distributions extend to x = 1.22 
corresponding to the r) threshold, or to s = 1.31, corresponding to the 27r threshold for D^''^ . 

to X = 1 — r. The size and shape of the smearing may affect the various integrals B^}'"^ differently and cannot 
be predicted, although we expect a small perturbative tail to emerge away from the endpoint. Clearly, smearing 
towards smaller Ef. tends to suppress the WA contributions to higher lepton energy moments, in which case we expect 
Bwa'^^''^ S Bwa'^- We also know that the WA contributions involve a SU{3) flavor singlet in the final state, and may 
be confined below the two-pion or rj thresholds which, for D decays and rric = 1.4GeV, correspond to a; < 1.31 and 
X < 1.22, respectively. For Dg decays the allowed region is slightly larger. Indeed, the WA distributions spread over 
a region of electron energies of 0{Aqcd)- In Fig- 2 we show four possible forms of the WA distribution as a function 
of X, assuming only positivity and the presence of a tail towards lower x. Two of the forms displayed stop at the rj 
threshold, the other two extend to the two-pion threshold. 

The endpoint behavior of WA can also be parameterized in terms of matrix elements of local four-quark operators 
of higher and higher dimension |19j . which determine the WA distribution. At leading order in the heavy quark 
expansion, however, the parameters -B^^*, i.e. the total moments of the WA distribution, are all determined by the 
dimension six term, and therefore they are all equal to B^^}^. The smearing is therefore a power and a^-suppressed 
effect in B^}'^, but it may be phenomenologically important despite the formal suppression. 

We will be primarily interested in a determination of the leading matrix element, B^}^ , i.e. the zeroth moment of 
the WA distribution, from b!^}'': the non-negligible smearing of WA may lead to a model-dependent dilution. In 
order to quantify this effect on the parameters -Bwl"^, we have considered a number of distributions like those in Fig. 2. 
In general the dilution is a moderate effect: b'^1'' > 0.8 B^^]'', > 0.7B^i^ and of course for WA distributions 

localized to the right of the partonic endpoint, B^^^^^'^ > i^wi'^- However, distributions characterized by a longer tail 
will lead to a stronger dilution, which may also be enhanced in the normalized moments. In the next Section, when 
we try to extract information on the WA from the leptonic moments, we will carefully take this effect into account. 

The WA contributions to the decays of charmed mesons involving different spectator quarks can be parametrically 
decomposed as follows: 



where 9c is the Cabibbo angle. In the isospin limit, the number of independent contributions is reduced by the 
identification B^}'^{D°) = B^l''^ {D~^). Any difference in the moments between the and the can only be due 
to isospin violation or Cabibbo suppressed contributions. Since the data do not show any significant difference in the 
rates or the moments of D° and D'^ in what follows we will neglect such Cabibbo suppressed effects. We are then left 
with just two distinct contributions, which can be identified with the valence and non- valence WA terms involving the 
s quark. The first one, i3wi'^(Z3s), only contributes to the Dg decays, while the second one, B^}'''{D), contributes 
equally to and decays. In the flavor SU{3) limit, the two contributions correspond to the sum of isotriplet 
and isosinglet, and to the isosinglet contribution, respectively. 




(10) 
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IV. RESULTS AND DISCUSSION 



Our first task will be to check whether the OPE at 0(l/m^) describes the experimental data in a satisfactory 
In the pole mass scheme with rric = 1.6GcV, r = 0.005, and /iwA = O.SGeV we have 



{E') = {E')o 



1 - 0.72 as - 0.29 a^^o - 0.60 fi% - 0.20 nl + 0.42 

.07/1^ + 0.20 - ^ 



1 - 0.03 a^- 0.03 a^/3o-0 
1 -0.07as - 
1 - 0.09 a, -0.05 



0.38 pis + 805^1 
as - 0.05 al(3o - 0.14 /i^ + 0.52 ^1 + 3.5 p% + 0.66 pig + 2045^, 



.05al/3o - O.Upl + 1.7 pI + 9.4 pf, + 1.4 pi^ + 641B^i 



way. 

(11) 
(12) 

(13) 
(14) 



where /3o = 9, Fq, {E)q, {E'^)q, {<^e)o are the tree-level results, and all coefficients are in GeV to the appropriate 
power. We immediately notice that the lepton moments receive smaller perturbative corrections in the pole mass 
scheme than the total rate, while they are very sensitive to WA and p^,. Since the uncertainty on the value of p^ to 
be employed in Eqs. (11 14) is sizable, it follows from the previous section that the choice of has a strong impact 
on the final error on WAT^Indeed, the simultaneous use of all three observables can reduce this ambiguity. As for 
the TOc dependence, since {E") scales like to", the lowest moments are less sensitive to the value of the charm mass 



than the width, which is proportional to to^ . We have also introduced effective WA parameters B^l ~ B 
B 



IB. 



(0) 



(2) 



B 



(2) 



2 p,(0) 



,(0) 

■^WA 



— 3.8i?wA 



o o d(2) 

O.ZJ^wA- 



ny 

WA — J-'WA ^J-'WA: J^WA — -Ov 

In order to avoid the well-known problems associated to the pole mass, semileptonic decays are usually described 
in terms of threshold heavy quark masses, like the kinetic mass, the PS, or the IS mass, see [Ij for a review. In the 
following we focus on the kinetic scheme |28| . For the charm mass and the non-perturbative OPE parameters we use 
as initial inputs the results of a global fit to i? — >■ Xcli' and radiative moments [29] , 



TOc(lGeV) = 1.16(5)GeV, 
pis(lGeV) = -0.19(8)GeV^ 



At^(lGeV) = 0.44(4)GeV^ 
p|5(lGeV) = 0.19(2)GeV^ 



The above values refer to expectation values of local operators in the B meson in the kinetic scheme with the IR 
cutoff scale pkin = IGeV. Up to power corrections, the D meson expectation values can be identified with the B 
meson ones, but pkin — IGeV is too high compared to the charm mass and we need to run the parameters to a lower 
IR scale above Aqcd- 

We evolve the OPE parameters including the charm mass down to pkin — 0.5GeV using 0{al) expressions. The 
expectation values Pq and p^g do not run with the kinetic scale. The running to such low scales induces significant un- 
certainties in the other parameters, of both perturbative and non-perturbative origin, which we estimate by varying the 
scale of as. We finally adopt TOc(0.5GeV) = 1.40(7)GeV, p^ (o.SGeV) = 0.26(6)GeV2, p|3(0.5GeV) = 0.05(4)GeV3. 
Of course, the charm mass determination we use is not the most precise. For instance, Ref. |30| reports a very precise 
value in the MS scheme which is consistent with |29j . However, any scheme translation would increase the error 
significantly. In the end, we do not expect an improvement even in the case of the width. In addition, we assume 
that the large uncertainty on p^g as extracted from the B fit dominates over the perturbative and non-perturbative 
corrections to its value. In the case of pg, on the other hand, we use a conservative approach with the central value of 
[5T] and a large error p^ = 0.35(10)GeV2. For the strange quark we use its MS definition, TO5(2GeV) = 0.105(2)GeV 
[32]. We evaluate the corrections to the total rate and the first two leptonic moments at pwA — 0.8GeV as 

^k^n = 1.2(3)10"^^Gev|l + 0.23a, +0.18a^/^o- 0.79 PG-0.26p^ + 1.45 p|,+0.56pi5 + 120B|!?i|, (15) 

iEe)k^n = 0.415(21)GeV { 1-^0.03 a, 0.02 af/^o- 0.09 p^ -I- 0.26p2 -i-2.7p|)+0.44pis-)- 203BiJi} , (16) 

0.192(20)GeV^ |l + 0.001 a, + 0.02a^jQ ~ 0.18 p^ + 0.68p^ + 6.6p% + 0.99pis + 307B^1| , (17) 

019(2)GeV^ |l - 0.53as - 0.17a^/3o - 0.18pG + 2.2p^ + 17p^ + 2.lpls + 961i?^2| , (18) 



iEe)kin 



(^E,k 



= 0.1 



where all coefficients are in GeV to the appropriate power and the tree-level results include the error due to their mass 
dependence. In the first and the second leptonic energy moment, the dominant source of uncertainty is by far p^,, 
while in the total rate the charm quark mass uncertainty dominates. In view of the accuracy relevant to our analysis, 
we can turn a blind eye to the bad convergence of the expansion for the width. Varying all the input parameters as 
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described above and as by ±20% around as(TOc) = 0.36 to account for higher order perturbative effects, we compare 



the theoretical 


expressions with eqs. 


(2,6) and we 


extract the following 


values of B^l{Dq) and B^^'^^Dq) 




s^i' (^°) = 


= -0.001(3)GeV^ , 




= -0.001(3)GeV^ , 


B^-^^'^^^iDs) 


= -0.002(3)GeV^ , 


(19) 




= -0.0001(6)GeV^ , 




= -0.0001(6)GeV3 , 


BwA^{Ds) 


= -0.0001(6)GeV^ , 


(20) 




= -0.0001(10)GeV^ , 




= -0.0002(10)GeV^ 


-Swa'*(-Ds) 


= -0.0001(10)GeV^ , 


(21) 




= -0.0000(7)GeV^ , 




= -0.0000(7)GeV3 , 


B^^l'^iOn 


= 0.0001(10)GeV^ 


(22) 



At the estimated theoretical precision, all the extracted WA contributions are consistent with zero. The obvious 
implication is that the OPE describes all the data reasonably well. We also repeated the exercise at different kinetic 
scales, finding in general similar results, despite the high sensitivity to accidental cancellations. For larger Hkin the 
size of perturbative corrections to the rate and their associated error increase rapidly as noticed in [34; . However, the 
moments are less affected and one obtains consistent WA estimates for kinetic scales as large as fikin — 0.8 GeV with 
similar errors. We have also repeated our analysis using the IS scheme for the charm mass |33j and fikin = for the 
other OPE parameters. The apparent convergence of the perturbative expansion improves for the rate and variance, 
but the estimates and errors of the WA contributions are very similar. 

Eqs.( 19|p2 1 show no departure from the SU{3) flavor symmetry, but the errors are strongly correlated among 



different mesons. On the other hand, in the ratios of Dg and semileptonic widths and leptonic moments, the 
dominant rric and all the perturbative corrections cancel, while only isotriplet WA and SU (3) breaking contributions 
to the OPE parameters remain. We estimate the latter conservatively as 20% of the central values and add the resulting 
errors linearly to account for possible correlations. Comparing these expressions to the experimental values we obtain 
AB^l'" = B^°l''{Ds) - B^^t'iD") = -0.0014(12) (5)GeV^ from the ratio of rates, AB^^t' = 0.0000(3) (l)GeV^ 
from the ratio of the first moments, and ABi^l'^ = 0.0000(4) (2)GeV'^ from the ratio of the second leptonic energy 
moments, where we show the theoretical and experimental uncertainty contributions in the first and the second 
brackets respectively. We observe a mild indication of a nonzero negative isotriplet iJwi'* contribution. This indication 
is however not really significant at present. 

So far wc have parameterized the WA contributions to the individual moments, and we have learned that the OPE 
gives a satisfactory description of the data even in the absence of WA, with the possible exception of AB^^Jf . As 
already mentioned in the previous S ection , we expect WA to be concentrated at the end-point, in which case we can 
combine the ZJ^'"*" results from eqs. (20 22 1 to obtain the isosing let B^l' = -0.0003(15) GeV^. In the worst case, the 



reduction in sensitivity due to the WA smearing functions shown in Fig. 2 is about 60%, which leads to 

B^^ = -0.0003(25)GeV^ , (23) 

where we have combined the dilution error linearly, without changing the central value. We regard this as sufficiently 
conservative. We can also extract AB^^I''* from AB^a^-*'*: in the absence of WA smearing they both give the same 
result, AB^t' = 0.0000(7) (3), while the worst dilution scenario gives AB^l' = 0.0000(12) (3), which is compatible 
with the value extracted from the ratio of widths. 

We can also look for indications of WA dilution by extracting B-^a^'^-* from a combined fit to the moments and the 

total rates (recall that B^^^ are linear combinations of -Bwa^'^"*)- The results for 1)° and Dg decays are shown in 
figure [3] and are consistent with B^l « B^l w -Bwl, but the errors are too large to draw a conclusion concerning WA 
dilution. 

Finally we study the correlation between and WA by considering a free parameter in the fit. In this way, 
the dominant source of uncertainty in the moments is removed. From the variance, first and second moments of 
wc find the following constraints (the results from Z?+ and Ds data are consistent) 

b'^^I'" + 0.017p|, = 0.0009(2)GeV^ , B^l'" + 0.021^1, = 0.0010(5)GeV3 , B^l' + O.OHpf, = 0.0006(4)GeV3 . 

Assuming no WA dilution one then obtains in each of these cases 

B^^ + 0.043p|) = 0.0021(5)GeV^ , B^^ + Ombpl, = 0.0016(9)GeV^ , B^^ + OmSp^, = 0.0014(10)GeV^ , 

and it is remarkable that they all agree quite well. This may be viewed as a mild indication that there is no 
significant WA dilution, or possibly that the dilution is similar in the three cases. Indeed, assuming vanishing WA at 
/iwA = O.SGeV all the moments can be reproduced by yo|3(0.5GeV) = 0.05(l)GeV^. 

In order to connect to B X^iv we refer to the WA matrix elements determined in charm decays as i?wA and 
consider their relation to those relevant in B semileptonic decays, i?wA, taking into account eq. \^\. In the heavy 
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— Xjl , I 

quark limit, /p ^ rrip |35j so that Bpi scale as constants with heavy quark mass, but recent lattice results give 
/£) 0.21 GeV, /b w 0.20 GeV [35]. We also neglect any evolution of the WA operators and write 

Sw'aI/^wa) = ^^B^a(Mwa) . (24) 

The parametric enhancement due to meson masses and decay constants is a significant factor of 2.5. Due to finite 
heavy quark masses, one might also expect additional power corrections, which spoil the exact scaling of i?wA between 
the D and B sectors. From eqs. ( [23p4| we derive a bound 



I^wa(Mwa = 0.8GeV)| < O.OOeCeV^ , (25) 

which holds for the non-valence contributions, although we stress that the data seem to prefer even smaller values. 
The valence contribution is more constrained from the ratio of the Dg and Z)*^ rates and leptonic moments, from 
which we expect 

- 0.004 GeV^ < ABl,^{fi^^ = O.SGeV) < 0.002 GeV^ (26) 

The bounds lead to a maximum 2% WA correction to the total rate of i? — > X^iv. In turn, this translates into an 
uncertainty of 1% on \Vub\ extracted from the total rate and from the most inclusive experimental analyses, like those 
that involve a lower cut on the invariant hadronic mass. Our bound on the WA expectation value improves on previous 
estimates. In [27], for instance, the maximum value allowed for i?^^(/iwA = IGcV) was as high as 0.020GeV'^, or 
-BwA(AiwA = O.SGeV) = O.OlSGeV^ 

Our analysis is compatible with the one of ref. [14 for the valence contribution, although we allow for a larger SU(3) 
violation. The analysis of ref. [T5] gives a determination of both valence and non- valence WA from the D semileptonic 
widths only. In our notation the results of ref. [15] correspond to a valence contribution AB^^ = — 0.0015(9)GeV^ 
and to an isosinglet contribution B-^^ — 0.0036(5)GeV^. The valence contribution is perfectly compatible with 
our estimate from the widths only. The non-valence contribution, on the other hand, has a tiny error and should be 
compared with B'^^ — — 0.001(3)GeV"^ that we extract from the widths only, see eq. (19). The difference is presumably 
due to various sources: i) the charm mass that we employ is higher than that employed in [15]: we start with the 
result of the global fit |29j . which is in excellent agreement with other rric determinations like charmonium sum rules, 
and evolve it to lower Hkin', H) our OPE parameters (taken from [29] and then evolved to lower /ifci„) are different 
from those of ref. [15], which may lead to sizable differences in charm but not in bottom decays; iii) the method of 
[15j implies the use of as(m,;,), and might underestimate the perturbative corrections. The discrepancy between the 
two determinations provides indeed additional motivation for using the moments to constrain WA, as they are less 
sensitive to power and perturbative corrections. In any case, we find remarkable that our determination from the 
moments is compatible with that from the widths and relatively stable wrt changes in the definition of the charm 
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mass or in fikin- Notice also that despite the larger value of B^j^ they find, the authors of [T5] end up with an estimate 



for the WA effects in the B — )■ X^ii^ rate which is very similar to ours, once the enhancement related to eq. (24) is 
taken into account. This is because our large uncertainty almost covers the discrepancy in the central values. 

V. SUMMARY 

We have performed an analysis of inclusive semileptonic D decays using the Heavy Quark Expansion. In addition 
to the total widths, we have used the Cleo data on the lepton energy spectra to compute the first few moments. 
The latter are quite sensitive probes of possible Weak Annihilation contributions, both in its isosinglet and isotriplet 
components, and determine very precisely a linear combination of the expectation values of the Darwin and WA 
operators. We have shown that the extraction of WA from the moments depends to the way WA is distributed in the 
lepton energy spectrum and we have taken this effect into account in our error estimates. 

Our analysis of Cleo data shows no evidence for Weak Annihilation, i.e. the OPE describes well the experimental 
results even in the absence of WA. There is a mild indication of non-zero valence WA in the ratio of the Dg and 
widths, but it does not seem to be supported by the ratio of the lepton moments in and 1?° decays. A possible 
explanation involves sizable SU{3) breaking in the matrix elements of the power suppressed operators and/or a WA 
contribution which is broadly distributed over the leptonic spectrum. 

We derive an upper limit on both valence and non-valence WA components, which allows us to put a bound of 2% 
on their relevance in the B — > Xuii' decay rate and even less for the isotriplet component. We look forward to the 
individual measurements of the B^ and B'^ charmless inclusive semileptonic decays, which could test our conclusions 
regarding the suppression of valence WA in these decays. 

Finally we note that the hadronic mass and moments in inclusive semileptonic decays of heavy quarks are quite 
sensitive to non-perturbative contributions, see e.g. |26j. Their measurement at Cleo-c or at the recently started 
BESIII experiment [37] might help us disentangle the various non-perturbative effects. 

Acknowledgments 

We thank Ikaros Bigi and Nikolai Uraltsev for enlightening discussions and various comments on the manuscript, 
the Galileo Galilei Institute for Theoretical Physics where this work has been finalized during the workshop Indirect 
searches for new physics at the time of LHC, and INFN for partial support during the completion of this work. JFK 
acknowledges useful discussions with Thomas Manuel and is supported in part by the European Commission RTN 
network, Contract No. MRTN-CT-2006-035482 (FLAVIAnet) and by the Slovenian Research Agency. PG is supported 
in part by a EU's Marie-Curie Research Training Network under contract MRTN-CT-2006-035505 (HEPTOOLS). 



[1] M. Bona et al. [UTfit Collaboration], arXiv:0908.3470' [hep-ph], see also http: //www.utf it . org/', for recent \ Vub\ averages 

see http: //www. slac . Stanford. edu/xorg/hf ag/ 
[2] M. Antonelli et al, arXiv:0907.5386 [hep-ph]. 

[3] J. Chay, H. Georgi, and B. G rinstein , Phys. Lett. B247 (1990) 399; I. I. Y. Bigi, N. G. Uraltsev, and A. I. Vainshtein, 
Phys. Lett. B293 (1992) 430 |arXiv:h ep-ph/9207"214| [Erratum, ibid. B 297 (1993) 477]; I. I. Y. Bigi, M. A. Shifman, 
N. G. Uraltsev, and A. I. Vainshtein, Phys. Rev. Lett. 71 (1993) 496 [arXiv:hep- ph/9304225]; A. V. Manohar and 
M. B. Wise, Phys. Rev. D49 (1994) 1310 |arXiv:hep-ph/9308246 . 

[4] I. I. Y. Bigi and N. G. Uraltsev, Nucl. Phys. B 423, 33 (1994) arXiv:hep-ph/9310285 . 

[5] R. D. Dikeman and N. G. Uraltsev, Nucl. Phys. B 509, 378 (1998) arXiv:he p-ph/9703437 ; I. I. Y. Bigi, R. D. Dikeman 

and N. Uraltsev, Eur. Phys. J. C 4, 453 (1998) |arXiv:hep-ph/9706520 . 
[6] M. S. Baek, J. Lee, C. Liu and H. S. Song Phys. Rev. D 57 (1998) 4091 |arXiv:hep-ph /9709386| ; H. Y. Cheng and 

K. C. Yang, Phys. Rev. D 59 (1999) 014011 [arXiv:hep-ph/9805222 . 
[7] M. Di Pierro and C. T. Sachrajda [UKQCD CoUaboration] , Nucl. Phys. B 534 (1998) 373 f arXiv:hep-lat/9805028| . 
[8] D. Becirevic, arXiv:hep-ph/0110124, 

[9] D. Pirjol and N. Urahsev, Phys. Rev. D 59 (1999) 034012 arXiv:hep-ph /9805488| . 
[10] N. Uraltsev, Int. J. Mod. Phys. A 14 (1999) 4641 arXiv:hep-ph/9905520'^ 



[11] I. I. Y. Bigi and N. G. Urahsev, Z. Phys. C 62, 623 (1994) |arXiv:hep-ph/93 11243 

[12] M. B. Voloshin, Phys. Lett. B 515, 74 (2001) arXiv:hep-ph/0106040 . ' 

[13] D. Becirevic, S. Fajfer and J. F. Kamenik, Phys. Lett. B 671, 66 (2009) lar Xiv:0804.Tf50| [hep-ph]]. 
[14] I. Bigi, T. Manuel, S. Turczyk and N. Urahsev, arXiv:0911.3322 [hep-ph]! 
[15] Z. Ligeti, M. Luke and A. V. Manohar, .arXiv:1003. 1351. [hep-ph]. 



10 



[16] D. M. Asner et al. [The CLEO Collaboration], arXiv:0912.4232 [hep-ex], 

[17] B. Blok, R. D. Dikeman and M. A. Shifman, Phys. Rev. D 51 (1995) 6167 [arXiv:hep-ph/9410293 



[18] A. Bazavov et al. [Fermilab Lattice and MILC Collaborations], PoS LAT2009 (2009) 249 arXiv :0912.522T] [hep-lat]] . 
[19] A. K. Leibovich, Z. Ligeti and M. B. Wise, Phys. Lett. B 539 (2002) 242 arXiv:hep-ph/0205148 . 

[20] A. Pak and A. Czarnecki, Phys. Rev. Lett. 100, 241807 (2008) arXiv:0803.0960 [hep-ph]] and K. Melnikov, Phys. Lett. B 

666, 336 (2008) arXiv:0803.0951 [hep-ph]]. 
[21] V. Aquila, P. Gambino, G. Ridolfi and N. Uraltsev, Nucl. Phys. B 719 (2005) 77 arXiv:hep-ph /0503083| . 
[22] M. Gremm and A. Kapustin, Phys. Rev. D 55 (1997) 6924 arXiv:hep-ph/9603448]. 



[23] B. M. Dassinger, T. Mannel and S. Turczyk, JHEP 0703, 087 (2007) arXivihep-ph/MTTTeS] . 

[24] P. Gambino and P. Giordano, Phys. Lett. B 669, 69 (2008) arXiv:0805.0271 [hep-ph]]. 

[25] A. Czarnecki and M. Jezabek, Nucl. Phys. B 427, 3 (1994) |arXiv:hep-ph/9402326 . 

[26] P. Gambino, G. Ossola and N. Uraltsev, JHEP 0509 (2005) 010 arXiv:hep-ph/050509T]. 

[27] P. Gambino, P. Giordano, G. Ossola and N. Uraltsev, JHEP 0710 (2007) 058 arXiv:0707.2493 [hep-ph]]. 

[28] I. I. Y. Bigi, M. A. Shifman, N. Uraltsev and A. I. Vainshtein, Phys. Rev. D 56 (1997) 4017 | ar Xiv: hep-ph /9704245| and 

Phys. Rev. D 52 (1995) 196 arXiv:hep-ph/94054lO . 
[29] O. BuchmuUer and H. Flacher, Phys. Rev. D 73 (2006) 073008 ' arXiv:hep-ph/0507253 , updated for the PDG 09 web 

update, http://www.slac.stanford.edu/xorg/hfag2/semi/winter09; P. Gambino and N. Uraltsev, Eur. Phys. J. C 34 

(2004) 181 | arXiv:hep-ph/0401063 . 
[30] L Allison et al. [HPQCD Collaboration], Phys. Rev. D 78 (2008) 054513 ar Xiv:0805.2999] [hep-lat]] . 
[31] N. Uraltsev, Phys. Lett. B 545 (2002) 337 arXiv:hep-ph/0111166 . 

[32] C. Amsler et al. [Particle Data Group], Phys. Lett. B667, 1 (2008) and 200 9 partial update for the 2010 edition. 
[33] A. H. Hoang, Z. Ligeti and A. V. Manohar, Phys. Rev. D 59, 074017 (1999) | arXiv:hep-ph/981 1239 | . 
[34] J. F. Kamenik, arXiv:0909.2755 [hep-ph]. 



[35] B. Grinstein, E. E. Jenkins, A. V. Manohar, M. J. Savage and M. B. Wise, Nucl. Phys. B 380, 369 (1992) [arXivihep- 
ph/9204207 . 

[36] V. Lubicz and C. Tarantino, Nuovo Cim. 123B, 674 (2008) [arXiv:0807.4605| [hep-lat]] . 



[37] D. M. Asner et al., arXiv:0809.1869| [hep-ex] 



